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Abstract 

In the normal phase (where no dynamical fermion mass generation occurs) of 
the D-dimensional quantum electrodynamics with Nf flavors of fermions, we derive 
an integral equation which should be satisfied by (the inverse of) the wave function 
renormalization of the fermion in the Landau gauge. For this we use the inverse 
Landau-Khalatnikov transformation connecting the nonlocal gauge with the Landau 
gauge. This leads to a similar equation for the running flavor number in the framework 
of the resumed Schwinger- Dyson equation. Solving the equation analytically 

and numerically, we study the infrared behavior and the critical exponent of the 
3-dimensional QED (QED3). This confirms that the flavor number in QED3 runs 
according to the (3 function which is consistent with the asymptotic freedom as that 
in 4-dimensional QCD. 



* To be published in Physics Letters B. 



Several years ago [0] running of the flavor number Nf has been investigated in the 
three-dimensional QED (QED3) based on the improved Schwinger-Dyson (SD) 
equation. It was shown that the parameter 1/iV/ exhibits the asymptotic freedom in 
QED3 as in non-Abelian gauge theory in 3+1 dimensions. This conclusion was drawn 
by solving the SD equation for the wave function renormalization for the fermion in 
the Landau gauge including the vertex correction. In order to arrive at the above 
conclusion, however, we needed to take several approximations. They are reexamined 
and some of them are criticised in p[. In the previous work P] the running of the 
flavor number was translated into the behaviour of the nonlocal gauge which can be 
derivable for a more general class of the vertex function than that considered in 
This work has resolved the criticism raised by pf. 

The purpose of this Letter is to rederive the running of the flavor number based on 
a different method which is superior to the previous one in the sense that the new 
method is free from the various ambiguities and most of the approximations taken in 
the previous paper . This result enables us to understand the normal phase of high- 
Tc superconductor [Q as a non-trivial IR fixed point, see p. The equation derived in 
this Letter provides a useful tool to study, among other things, the existence of non- 
trivial quasi infrared fixed point and slowing down of the running flavor number in 
the intermediate region, which is responsible for the deviation from the Fermi-liquid 
behavior in the normal phase of high-Tc superconductivity. 

We focus on the SD equation for the full fermion propagator, S{p) = [A{p'^)p — 
in D-dimensional QED. The SD equation for S{p) is decomposed into a 
pair of two coupled nonlinear integral equations for A(p^) and B{p'^). In general, this 
procedure leads to rather complicated coupled equations. A method of avoiding the 
coupled equation is to find a situation such that the wave function renormalization 
disappears, i.e. Aip^) = 1 (no wave function renormalization). This is realized in 
the framework of SD equation. Q Actually, under the bare vertex approximation, the 
fermion wave function renormalization becomes 1 identically in the Landau gauge in 
any dimension D > 2, if the bare photon propagator is adopted in the SD equation of 
QED (i.e. in the quenched approximation), see e.g. Ref. This extremely simplifies 
the actual analysis of SD equation, since we have only to solve a single equation for 
B (See Ref. ^ for more advantages of the nonlocal gauge). This scenario has been 
extended beyond the quenched approximation 0, ||, |^, |T0|, |[ . Especially, it has been 
shown 1^ that in D-dimensional QED we can choose a covariant nonlocal gauge in 
which A{p'^) = 1 follows as a solution of the coupled SD equations, if the vertex 
function belongs to the following class: 

T,{p,q;k)=^,G{p',q',k'), (1) 

where G is an arbitrary function of the arguments except that G does not depend 
explicitly on /c^ (Note that implicit dependence on fc^ through A{k) is allowed [Q). 
Such a nonlocal gauge [| ^ = ^{k) (i.e. momentum-dependent covariant gauge fixing, 

^ In perturbation theory this is possible order by order in the couphng constant. 
^ This is obtained by integrating the first order differential equation satisfied by ^. Therefore, we 
need a boundary condition, see Ref. 0, ||, ||]. 
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see Ref. ||T0| , ^ for the Lagrangian formulation) is given by 



m :^ mil - n(kye] ^ 1 + /; (2) 

when the full photon (or gauge boson) propagator has the form 
D^^[k) = Drik ) U^^. — + 



k"^ J k"^ k"^ ^ 

«-(*^) - W^Y 

where T^{k) is the photon vacuum polarization function. In fact, in the quenched 
(n(fc) = 0) and ladder {G = 1) approximation, eq. (@) recovers the Landau gauge 
^{k) = for D > 2. Therefore, eq. (Q) is indeed a generalization of the quenched 
case. As a result, we have only to solve a single equation for B: 

B{P) = ^o + e J j^^^-^^^^DAk)[D-l+ak)]. (4) 

The mass function M(p) := B{p'^)/A{p'^) reduces simply to B{p^) in the nonlocal 
gauge (|^). The SD equation (without IR cutoff) has been solved for QED3 ana- 
lytically and numerically in |p, p!0| . 

A question is how one can recover the Landau gauge result in this scheme. As 
remarked in the previous work 0], the inverse Landau-Khalatnikov (LK) transforma- 
tion Dill can recover the full fermion propagator S{x), the full gauge boson propagator 
Dfj,u{x) and the full vertex function V^(a;, z) in the Landau gauge from those in the 
nonlocal gauge (H). 

The LK and the inverse LK transformations leave the system of the SD equa- 
tion form invariant. Usually the LK transformation is used to transform the results 
obtained in the Landau gauge into those in other gauges. This is because, in many 
cases, the calculation in the Landau gauge is easier than other gauges. (This is not 
necessarily true in the study of SD equation.) However, it is easy to see that the LK 
transformation is invertible, see Ref. [|12[. Therefore the inverse LK transformation 



enables us to transform the results in the arbitrary gauge (even the nonlocal gauge) 
into those in the Landau gauge. 

For example, in the D = d + 1 dimensional QED the fermion propagator Snlg{x) 
in the nonlocal gauge is connected with that Sl{p) in the Landau gauge (in configu- 
ration space) 0. 

SLix) = e-^^^^SNLGix), (5) 

where 

A(x):=e2/^(e^'=--l)/(fc), f{k):=^-^. (6) 
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Especially, in the normal phase (with no fermion mass generation, i.e. B{p^) = 0), 
Snlg{x) agrees with the bare fermion propagator S'o(x): 



^'^""^ = I (0^^"'^ = nVoix), Voix) = (7) 
Therefore, the momentum-space propagator Sl{p) in the Landau gauge 

Sl{p) = -— ^ (8) 

is obtained by applying the Fourier transform to (|^). Thus the wave function renor- 
malization in the Landau gauge is obtained in the normal phase |]^, 

Z{p) ■= A~\p^) = I J rf^a;e^f-^e-^(^)(p ■ x)Vo{x). (9) 

However, in order to obtain the wave function renormalization in the Landau gauge 
according to (H), we need to perform the Fourier integration twice (in (|]) and (P)). 
Although this is not impossible in principle, it becomes rather laborious to perform 
it when the nonlocal gauge ^{k) is a complicated function of k (This is indeed the 
case of QED3 in question). 

In this paper we propose another method which is superior to the above method 
(^ as shown in what follows. A basic observation is that in the normal phase the 
function Z{p) := l/Ai^p"^) in the Landau gauge satisfies the following integral equation, 

Z{P) = l + ^e' [ fl''^\~'^ fiP - Q)Ziq), (10) 



(27r)^ q 

where f{k) is given by (^. This equation is derived as follows. Note that in the 
normal phase Snlg^x) is equal to the massless bare fermion propagator 5*0(0;) which 
satisfies WSq^x) = 6^{x) where d := •y'^d^. Operating the differential operator d to 
both sides of (H), therefore, we obtain 

^^SL{x) = 5^(x) - ze'SM J ^kf{k)e-^'-^, (11) 

where we have used a fact A(0) = 0. Moving from the configuration space to the 
momentum space, we obtain 

pSl{p) = l + ze'J ^^Sl{p - k)kf{k), (12) 

since the Fourier transformation changes the simple product in coordinate space into 
the convolution in momentum space. By substituting (|^) into ([l^) , we arrive at the 
result eq. ([T0|) after the change of variable q = p — k. 

We can interpret the integral equation ([10|) as a SD equation derived from the 
fermion propagator S itself, not S~^. The conventional SD equation derived from 
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S^^{p) = A{p'^)p — B{p^) in momentum space leads to the SD equations for A and 
B, which are nonhnear integral equation in A and B. The exact solubility of eq. ([TO|) 
is due to the linearity of the equation in Z{p). It is worth remarking that the SD 
equation (|1^) can be derived from another type of the vertex function |jl2[. See 
Ref. [T^ for more details. 

In order to solve the equation (0), then, we move to the Euclidean space and 
separate the integration variables. Then we arrive at the equation, 



rA 

Z{p) = l-e^j dqZ{q)LD{p,q), 



(13) 



where p := \/]?, q := \f(^ and the kernel for D > 2 is given by 



D-3 



^ r d^sin^-^^ipqcos^ - q^)f{^J{p - qf). (14) 



Here we have introduced the ultraviolet (UV) cutoff A and the infrared (IR) cutoff e. 

In the following we restrict our arguments to the multiflavor QED3. In QED3 
the dimensionful quantity a can be used as an effective UV cutoff of the momentum 
integral, which is supported by the Wilsonian approach of the renormalization group, 
see and references therein. Hence we put A = a hereafter. (Application of eq. ([I0|) 
to other models will be presented elsewhere [0.) In the normal phase, the vacuum 
polarization (at one loop) reads 



-ak. 



a 



(15) 



Note that this quantity is gauge invariant. 

The nonlocal gauge in the absence of IR cutoff 



Then we obtain 



fik) 



k"^ + ak 
'~k^ 
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k'' - 2ak + 2a^ln 1 + 



is given by 

k^ 



a 



2a 2^2 , 



1 + 



k^ 



a 



k^ + ak^ 

This should be substituted into the kernel given by 

L3{p,q) '■= ^ d^sm^{pq cos^ - q^)f{\l{p - qf). 



(16) 



(17) 



(18) 



Note that the kernel L3{p,q) is not symmetric in p,q, i.e. Ls{p,q) 7^ L^{q,p). The 
angular integration in the kernel L^{p., q) can be done exactly, see Appendix. However, 
the expression is rather involved. The full analysis of equation (|l^) with an exact 
kernel L^{p, q) will be done using the numerical method later. 

First, we give an analytical solution by taking into account the leading term only 
in the IR expansion: 



f{k) 



E2(-i) 

n=0 



1 



, n + l _ 2 
n + 3 3a/c^ a^k"^ 



(19) 
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and then discuss the effect of the other terms. It should be remarked that the factor 
2/3 in the leading term of f{k) comes from the IR value ^(0) = 2/3 of the nonlocal 
gauge ( ]TB| ) in the absence of IR cutoff ^ : 



Up to the leading term in (|1^), the kernel reduces to a very simple form: 



1 1 

TT^ Sag 



9{q-p), 



(20) 



(21) 



where 6{x) is the Heaviside step function. The integral equation (|13]) with the kernel 
(PI]) can be converted to the boundary value problem of the first order differential 
equation: 



Z{p) dp 



-7, 7 



37r2A^ 



/ 



(22) 



The general solution of the differential equation is given by Z{p) = Cp^'^ with a 
constant C. Imposing the UV boundary condition Z{a) = 1 at p = a, the solution is 
uniquely determined: 



Z{p) 



7 



37r2iV 



/ 



Note that Z{j)) ] +oo (or A{p^) [ 0) as p J, 0. The solution 
directly using the iterated kernel as follows. 



Z{p) 



°° 1 a 
1 + E4(7ln-^ 

n=l^' P' 



exp 



7 In 



a 



P 



(23) 

3D can also be obtained 
(24) 



which agrees with (P5|). The wave function renormalization in the IR regime ( ^5] ) has 
been conjectured based on either theoretical speculations in |T^, |15|, or in numerical 
treatment Jl^. Quite recently, essentially the same result as ( p3|) was analytically 
obtained independently in ||T^ by calculating (H) following the suggestion 0. ^ 

The validity of the IR expansion used above has been confirmed in the supercon- 
ducting phase where dynamical mass generation occurs [10[- However, the validity 
in the normal phase has not been confirmed yet. Note that the solution (^) is inde- 
pendent of the IR cutoff e due to the special nature of the kernel (^). This property 

^ Using the kernel K{p,q) = —e'^L''^\p,q) = i"//q)0{q — p), Eq. ( p^ ) is replaced by Z{p) = 
1 + dqK{p, q)Z{q). Then the solution is given by Z{p) = 1 + ' where {Kn ■ l)(p) = 

dqK(p,q){Kn-i ■ !)(<?) for (n > 2), and {Ki ■ = dqK{p,q)l{q) = 7ln(a/p). In this case, 
(if„.l)(p) = ^(lnA)n. 

^ Note that the prefactor in Z{p) is exactly 1 irrespective of Nf in (p^), while the complicated 
dependence on Nf appeared in p8|. 
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is not preserved if we include the next-to- leading term in (p!9D, since the kernel is 
modified as 



An' 



' -4 

3aq 



1 — 

2pq 



In 



p + q 
\p-q\ 



(25) 



Therefore, the analysis at least up to the next-to-leading order is necessary to see a 
possible effect of the IR cutoff, in light of the discussion on the cutoff dependence 
given in where the ratio a/e plays an essential role in the critical behavior. 

An advantage of the method proposed in this paper is that we can take into 
account all the terms up to any order in a very similar way to the conventional SD 
equation. This will be indispensable for studying the intermediate region: e <^ p <^ a, 
which is responsible for the deviation from the Landau Fermi-liquid behavior and the 
existence of the non-trivial quasi IR fixed point [^]. 

A further advantage of this method is that we can directly write down the equation 
for the running flavor number (as a coupling constant) which was introduced in 



9Rip) ■= ^oM(/) = goZ{p), go := (26) 

where go is the bare coupling constant. Eq. (|13|) implies that the running flavor 
number in QED3 obeys 

gnip) = go- -jgo / dql^^p, q)gR{q). (27) 
4 JO 

Up to the leading term, the running coupling is obtained as 

p\-9o/^ / go p 



gnip) = go[-j = ^oexp (^-y In- ) . (28) 



The study of the equation ( P7| ) in the intermediate region will shed light on the 
existence of the non-trivial quasi IR flxed point of QED3 through the /3-function 

Pign) = j^gnit), t:=ln(p/«), (29) 

as well as the slowing down of the running in the intermediate region. 

In contrast to the optimistic observation of [113] , it is by no means obvious whether 



or not the leading results survive the inclusion of the higher order terms. Therefore, 
we have solved Eq. (|13|) for D = 3 numerically for the three kernels: the leading 
(pTD, the next-to- leading ( pSj) and the full kernel (|35D as given in the Appendix. The 
numerical result is given in Fig.l when Nf = 5 which lies in the normal phase (see 
the flnal comment). Fig.l shows that the solution (|23| ) in the leading order is a rather 
good approximation of the full solution of Eq. (0) in the IR region p/a < 0.1. In the 
IR region the full solution is approximated by (^) with the same 7. The IR critical 
exponent 7 deflned by (^2]) is nothing but the anomalous dimension. Therefore it 
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turns out that the result of |T8| is in fact correct with regard to the critical exponent 
7, the result for A being changed merely by a multiplicative constant. 

In conclusion our method confirms that the leading solution (^) for the running 
flavor is also a good approximation of the full solution of (0). We find that Eq.(]28|) 
is rewritten as 

gnip) = 9o[l - ^ In ^ + 0{g',)] = go/il + f In ^) + Oigl), (30) 
which implies the /3-function: 

f3{9R) = -l9l + 0{gl). (31) 

Thus QED3 exhibits the asymptotic freedom for the coupling go, at least for small 
go, as claimed in 

In the nonlocal gauge, the running coupling is related to the nonlocal gauge as H] 



NLGf,\ 

9r W = 90 



1 + liik') 



t := ln(A;/a). (32) 



This means that the growth of the running coupling is cut off at the IR limit k ^ 0, 
since ^(0) is finite irrespective of the existence or absence of IR cutoff, in sharp contrast 



to the result (|28|), i.e. gnip) t +00 as p | 0. Therefore, it is interesting to see the effect 
of the covariant IR cutoff in the Landau gauge and to see whether the situation 
change essentially or not. In a subsequent paper |T3| we will give more detailed 



analysis on the running coupling and the fixed point in the presence of covariant IR 
cutoff. For this we need to extend the LK transformation so as to include the massive 
gauge field. 

Finally, we comment on the superconducting (chiral- symmetry breaking) phase 
where the dynamical fermion mass generation occurs. It is worth remarking that 
Appelquist et al. [|19] assumed A{p'^) = 1 in the Landau gauge from the beginning 



based on the naive argument. This became an origin of the criticism later by 

Pennington et al. [0. Moreover it has been suggested that the discrepancy may 
be explained through the effect of the IR cutoff However the analyses H, 
based on the nonlocal gauge have confirmed most of the original results fl^. Among 
other things, they proved that there exists a finite critical number of flavor = 
128/(37r^) = 4.32 below which {Nf < Nj) the chiral symmetry is spontaneously 
broken and the fermion mass is dynamically generated, and that the scaling law is 
given by the essential singularity type (at the critical point Nj) for the dynamical 
fermion mass and the chiral order parameter The resulting value of A'^^ agrees 

with that of the improved calculation incorporating 1/A^| corrections by Nash ||20|| . 
Note that the above result in the neighbourhood of the critical point does not depend 
on the explicit choice of G, as long as the vertex is consistent with the Ward-Takahashi 
(WT) identity, see §. 
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A Kernel 



By changing the integration variable from ^ to r := \fW = i/p^ + q'^ — 2pq cos 6, the 
kernel (111) can be rewritten as 



1 rp+q 

4vr^^3(]9, q) = — dr{p^ - q^ - r'^)rf{r). 

Zpq J\p-q\ 

Integrating the RHS term by term using the expansion (|19D of f{k), we obtain 



(33) 



^n^Up,q) = — 



2E 

n=0 



-1)" (p + g)"+^ -\p-q 



\n+l 



ra + 3 



a 



n+l 



p + q 



l)"(n + 1) {p + qy-^ -\p-q 



|n-l 



a 



71+1 



p- - q- 



a-" 



In 



\p-q\ 



(34) 



Here note that the last logarithmic term corresponds to the n = 1 case which is 
ommited in the summation in the second term on the RHS of the above equation. 
Up to = or n = 1, this reproduces (^TJ) or (P^D, respectively. 

On the other hand, using (|l^ for /(fc), we arrive at the complete kernel: 



47r'L3(p,g) 



a{p^ — q^) 
2pq 

ipq 



1 



\p - q\^ (p + q)^ 
1 1 



_ p2 _^ g2 



2apq 



\p-q\ (p + q) 



+^ln 



2a^pq 



a + p + q \p — q\ 
a+\p - q\ p + q 



a 



2pq 



p + 3g, a + p + q P - 3q , ct + \p - q\ 
In -T7 m 



{p + qf 



a 



{p-qf 



a 



(35) 



Indeed, this agrees with (p^ after expanding (^) in powers of 1/a. 

This should be compared with the kernel in the SD equation @) for B in the 
nonlocal gauge. 



g{k) ■.= DT{e)[2 + l{k' 



k"^ + ak k"^ 



' 2a / k^ 



(36) 



This leads after angular integration to the kernel |10 

a p + q — \p — q\ 1 



M'K;{p,q) 



In 



a + p + q 



pq\p-q\{p + q) pq a+\p-q 



+- 



a 



pq 



{p + q)'- 



In 



a + p + q 



a 



\p-q\ 



In 



a+\p- q\ 



a 



(37) 
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If we use the expansion 

m = E 4(-ir = A - 4 + 0{k), (35 

the kernel (^) reads 

1 4 ~ (-l)"(n + 2) + - 



8vr^ir3(p, g) = —rjdrr,(r) = - ^ 



p + g— |p — g| 9 1 



3a 



2j9g 8 a 



+ 0(1K). (39) 



It has been analytically shown that the inclusion of the next-to- leading term in 
the kernel does not change the critical behavior near the critical number of flavors. 
Subsequently, this has been confirmed by the numerical study using the complete 
expression (0) in Ref. |TD|. 
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Figure Caption 

Fig.l: Wave function renormalization Z{p) as a function of p/a for the leading 
1]), the next-to-leading (p5|) and the full kernels ([35| ) when Nf = 5. 
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